Dijet induced collective modes in an anisotropic quark-gluon plasma 
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We discuss the collective modes due to the propagation of two oppositely moving relativistic 
jets (dijet) in an anisotropic quark-gluon plasma(AQGP) and compare the results with the case of 
single jet propagation. For the sake of simplicity, assuming a tsunami-like initial jet distribution, we 
observe that the dispersion relations for both the stable and unstable modes are altered significantly 
due to the passage of dijet in comparison with the case of single jet propagation. It has been 
further demonstrated that the growth rate of instability, due to introduction of dijet in the system, 
increases compared to the case of single jet case. As in the case of single jet propagation, the 
instability always grows when the jet velocity is perpendicular to the wave vector. We, thus, argue 
that the introduction of dijet in the AQGP, in general, leads to faster isotropization than single jet 
propagation. 
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I. INTRODUCTION 

One of the goals of the ultra-relativistic heavy-ion col- 
lision experiments at BNL RHIC and at CERN LHC is 
to produce a deconfined state of QCD matter commonly 
known as quark gluon plasma (QGP). The properties of 
such a system, if formed, can be studied through various 
probes such as electromagnetic probe, jet quenching etc.. 
According to QCD thermodynamics studied on lattice, 
this novel state of matter is expected to be formed at 
temperatures of the order of T ~ 170 - 200McV. The jet 
of particles produced initially due to hard collision, when 
passes through such a system loses energy and results in 
decrease of high energy hadrons produced due to the jet 
fragmentation. This phenomenon is commonly known as 
jet quenching. Moreover, the passage of jets influences 
the collective modes of the system [lHa] ■ 

The effect of single jet propagation on the plasma col- 
lective modes has been studied in Refs [U-IH where it has 
been assumed that the plasma is isotropic in momen- 
tum space and the jet is described by a tsunami-like jet 
distribution. In a very recent work the present authors 
have relaxed the condition of isotropic QGP and studied 
the characteristics of both stable and unstable collective 
modes of AQGP due to the passage of relativistic jet us- 
ing transport theory Q . It is shown Q that the disper- 
sion relations are modified quite substantially in compar- 
ison to the no jet case. On the other hand, the authors 
of Ref Q , have investigated the dijet induced properties 
of the unstable modes in isotropic QGP. Also the energy 
loss due to stream instabilities induced by dijet has been 
discussed using both tsunami-like and Gaussian distribu- 
tions for the relativistic jets within the same approach as 
in Refs. [IH3 . We extend the work of Ref. 0, [| to AQGP 
and study the characteristics of both stable and unstable 
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modes due to the passage of two jets moving in oppo- 
site direction. It should be noted that non-equilibrium 
jet of particles while travelling through QGP disturbs 
the plasma exciting chromomagnetic and chromoelectric 
modes of which some of modes might even be unstable. 
The most important among those are the modes which 
grow exponentially in time. Plasma instabilities due to 
jet propagation in AQGP have been the subject of ex- 
tensive studies as when the instability occurs, the kinetic 
energy of the particles is converted to the field energy 
leading to faster isotropization 6-9] of the QGP. This 
type of phenomenon occurs if it is assumed that the the 
hadronization time is greater than the time to generate 
the growth of the gauge fields. For the complete devel- 
opment of plasma instabilities, the time scale should be 
of the order of t ~ (6.7 — 12.5) /uJt, where uj t is the total 
angular frequency of the whole system 0, El] . 

We use the same simplifying assumptions as in Ref. @- 
IH to study the stable and unstable modes of AQGP 
induced by two counter propagating jets and compare 
it with that obtained in case of single jet [J] and no 
jet case [HI El- AQGP is realized in the very early 
stages of heavy-ion collision due to rapid longitudinal 
expansion. As a consequence, the cooling is faster in 
the longitudinal direction which results in (p 2 L ) « (pj). 
Such momentum-space anisotropy modifies the collective 
modes which have different behavior in comparison to 
that in isotropic QGP (see Refs. [HI [HI for details). 
Due to the passage of relativistic jets in a non-equilibrium 
plasma, with an anisotropic distribution in momentum 
space, the behavior of collective modes change as will be 
demonstrated in the following. 

The essential ingredient to study the collective modes 
in the composite system (AQGP+jets) is the transport 
equation which reduces to Vlasov equations for the time 
scale shorter than the mean free path time 0, HH • In 
our calculation we assume (i) weak coupling regime i.e. 
g « 1 and (ii) the interactions between jet and plasma 
is only mediated by mean gauge fields. The interactions 
between the jet particles and the plasma particles leads 



kinetic instability initiated either by charge or current 
fluctuations which lead to electric and magnetic instabil- 
ities. The latter is also known as Weibel instability fl3| . 
It has been demonstrated in Ref. [To|, [Hi that the growth 
rate of magnetic instability is maximum in the direction 
of anisotropy in momentum-space anisotropic plasma. So 
while studying the characteristics of the unstable modes 
we assume that the momentum of the collective mode is 
in the direction of the anisotropy. 

The organization of the paper is as follows. In section 
2 we mention the required formalism very briefly as detail 
has already been presented in [H-Q. In section 3 results 
for both stable and unstable modes will be presented fol- 
lowed by summary in section 4. 



II. FORMALISM 

To study the jet induced collective modes in AQGP 
we dwell on the transport theory described in detail in 
Ref. 0-0] ■ For the sake of simplicity, we assume two 
jets propagating in opposite direction described by the 
following tsunami-like distribution Q 



f jet (p)=nu°S^(p-Au). 



(1) 



Here n is related to the density of the plasma and u M 
is the four-velocity. The parameter A fixes the scale of 
energy of the particles. More realistic distribution, such 
as, a Gaussian distribution can be used to simulate the 
same effect 

Using Vlasov approximation one can write down the 
polarization tensor for particles species a as 0, 0, [HI [IH : 



rc(fc) 
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, a/q(p) 



p 



p.k 



(2) 



where a specify the quarks, antiquarks, gluons or par- 
tials of jet. The polarization tensor is symmetric, i. e. 
IP" (A) = IT^(fc), and transverse, fc^IP 1 " = 0. 

We assume that the ansatz for the anisotropic momen- 
tum distribution is given by fiol [ll| : 



(3) 



Here fi SO is the distribution function in isotropic case. 
iV(£) is the normalization constant which is equal to 
\J\ + £, n is the unit vector along the direction of 
anisotropy. The anisotropy parameter £ lies in the range 
— 1 < £ < oo. Using the above ansatz the spacelike com- 
ponent of the self-energy tensor can be written as (from 
Eq.©): 
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v j k l 
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is the isotropic Debye mass. Since the self energy depends 
on the four-momentum^^ 1 ) and the the anisotropic vec- 
tor (n^ = (l,n)), it can be decomposed in a suitable 
tensorial form in a covariant gauge jloL fl6j : 



n^'(fc) = a A 13 + /3B ij + 7 C' y ' + SD ij 



where 



A 11 



S ij - fcV/k 2 , 
k l k j /k 2 7 



= k l 



(6) 



(7) 
(8) 
(9) 
(10) 



where n l = A l ^v? which obeys n.k = and the structure 
functions a, /3, 7 and S can be obtained by appropriate 
contractions (see [1(1 EH for details). In the limit of 
£ — > 0, (isotropic case) 7 and S vanish whereas, a and 
(3 are directly related to the transverse and longitudinal 
components of the polarization tensor of the plasma re- 
spectively. 

In a similar way, the polarization tensors induced by 
the jets (moving in opposite direction) with the tsunami- 
like momentum distribution can be calculated to obtain 
the following expressions Q: 



(c 2 -k 2 V et ^ 



J''l 



k.v 



jet 



(U - k.V Jet ) 2 



u) + k.v jet (w + k.v jet ) 2 



where Vj et is the velocity of jet and w 2 et = ^ is the 
plasma frequency of the jet. Moreover, the polarization 
tensors due to the jets can be decomposed in the following 
way: 



ir/ et 2 = o"a« + p" 



(13) 



(14) 



where 



-[2 



L0 2 -k 2 



(UJ - kVjet COS jet) 
,2 „„ 2 , 



: v 2 sin 2 6 jet] 



2 2 V COS 6 j et — 1 
3et (u) - kv jet COS^et) 2 



(15) 



and 
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For numerical convenience, we introduce another param- 



eter r\ defined by r\ = w 2 et /w 2 where 



In the analysis of the collective modes of the composite 
system we arc interested in very short time scales where 
the collisionless approximation is justified. The effect of 
jet of particle is to induce the color fluctuations, which 
provide a contribution to the polarization tensor of the 
system. According to the linear response theory, the total 
polarization tensor of the whole system is given by the 
sum of polarization tensors of the plasma and the two 
jets. Thus we have, 

nr (fc) = n^(fc) + n^(fe) + n% t2 (k) (17) 

The dispersion relation of the collective modes of the to- 
tal system can be determined by solving the equation 

det[(fc 2 - lj 2 )S' 13 - k l v + rr/(fc)] = det[A- x (fc)p = o. 

(18) 

where, in the temporal axial gauge, the effective propa- 
gator A y is given by, 



A ij (k) 



1 



[(k 2 — uj 2 )5^ — fc'/V + nf (k)] 



(19) 



Therefore, using Eqs.(6), (11), (12), and (17) the effective 
propagator of the composite system can be written as 



A(fc) = A A A + (k 2 - uj 2 + a 



7 )A G B 



+ [(0 + 0' + 0" -u 2 )A G - A A ]C-<5A G E(20) 



with 



A7 1 = k 2 - uj 2 
Ac 1 



(k 2 - u 2 + a + a' + a" + 7 )(/3 + 0' + 0" 

kW. 



(21) 
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FIG. 1: (Color online) The dispersion relation for the stable 
A-mode for an anisotropic plasma with jet for different values 
of the anisotropy parameter £ = {1, 10}, r\ = 0.2, 8j et = 0, 
Vjet = 0.7 and 8 n = 
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III. RESULTS 



A. Stable modes 

The collective modes of the system comprising of the 
plasma and the jets are obtained by finding the zeros 
of Eqs. ((21]) and <(22j). For real values of u; > |k|, the 
dispersion relation for A-modes can be determined by 
finding the solutions of the equation: 

u)\ = k 2 + a{u> A ) + a'(uj A ) + a"(uj A ) (23) 

In case of G-modes we factorize A^ 1 as 

Ag 1 ^^-^ 2 -^) (24) 

where 

w%± = \{Co 2 ± \]W +4k 2 n 2 <5 2 ), (25) 



FIG. 2: (Color online) The dispersion relation for the stable 
A-mode for the system composed of anisotropic plasma with 
single jet and dijet for £ = 10, Vj e t = 0.7, 6 n — 0, 77 = 0.2 and 
9 jet ='{0, tt/3}. 

and 

Q 2 = a + a' + a" + + p + /3" + 7 + k 2 , 
il = a + a' + a" -0 - 0' - 0" + 7 + k 2 . (26) 

For real u> > |k|, the square root of Eq. ([23]) is always 
positive leading to two stable modes for G-modes. 

We vary various parameters introduced earlier to find 
the solutions numerically. We, in FigfTJ plot the results 
for the stable A-mode with single jet and dijet for various 
parameters explained in the figure. For completeness we 
have also plotted the results when there is no jet. We see 
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FIG. 3: (Color online) Same as Fig[2] for different values of 
6 n = {tt/4, 7r/2} and £ = 10, 9 jet = 0, f] = 0.2 and v jet = 0.7. 
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FIG. 4: (Color online) Same as Fig[2]for £ = 10, 6 n = tt/4 
0jct = tt/3, ry = 0.2 and « Je t = {0.5, 0.7} 



that the collective modes with the dijet differs reasonably 
from that with single jet and also with no jet. In order 
to see the the dependence of the dispersion relation on 
the angle of propagation of the jet with the wave vector 
the results for the A-modes have been displayed in FigJ3] 
for £ = 10, Vj e t = 0.7, 9 n = and rj = 0.2 where one 
can clearly see that the dispersion relation is sensitive 
to 9 j e t . Next we consider the variation of the dispersion 
relations for A- modes with the anisotropy direction 9 n . 
For fixed 9j e t the results are in shown in Fig(3J The dis- 
persion relation for dijet (jet) is not at all sensitive at 
low momentum and is marginally sensitive at higher mo- 
mentum. The important thing is to notice that the dijet 
results for fixed 9 n substantially differ from that with 
the single jet. The sensitivity of the dispersion relation 



FIG. 5: (Color online) The dispersion relation for the stable 
G+-mode for an anisotropic plasma with jet for 9„ — n/4, 
djet — 0, rj — 0.2 and Vj e t = 0.7. The left(right) panel corre- 
sponds to anisotropy parameter £ = 1(10). 
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FIG. 6: (Color online) The dispersion relation for the stable 
G+-mode for the system composed of anisotropic plasma with 
single jet and dijet for £ = 10, Vj e t = 0.7, 6 n = 7r/4, 77 = 0.2 
and different 8j et = {0, tt/2}. 



for the stable A-modes on the jet velocity is displayed 
in Fig |4] for 8 n = 7r/4. We see that the results changes 
marginally by varying the jet velocity both for single jet 
and dijet propagation. 

It can be checked that 7 and n 2 vanish identically when 
k||n. This leads to similar dispersion relations for G+- 
modes for 9 n = 0. First of all, we display the results 
for G+-modes in FigJS] for two values of the anisotropic 
parameter £. The results are quite sensitive to £ for fixed 
9 n . It is also seen that for fixed £ and 9 n the dijet re- 
sults differ substantially from the single jet case. We 
also find marginal dependence of G+-modcs on 8j et at 
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FIG. 7: (Color online) Same as FigH for£ = 10, v jet = 0.7, 
0jet = and t? = 0.2 and 9„ = {tt/4, tt/2}. 
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FIG. 8: (Color online) Same as Fig© for 6 n = tt/4, £ = 10, 
r? = 0.2, e jet = vr/3 and v jet = {0.7, 0.5}. 



very low momentum(see Fig|B]). In FigO the sensitiv- 
ity of the dispersion relation of the G+-modes with the 
anisotropy direction (9 n ) has been displayed. It is found 
that the stable G+-modes are quite sensitive to 9 n . It 
is also observed that the dijet results are quite different 
from the single jet case which is an important finding in 
the present work. We have also found that as in the case 
of single jet, dijet results are not that sensitive to the jet 
velocity (see Fig(5|). Finally, for stable modes, we con- 
sider the G_ -modes which have been shown in FigJ5]for 
£ = {1, 10}, 77 = 0.2, e jet = 0, v jet = 0.7 and 6 n = tt/4 
with and without jets. We find significant difference in 
the collective modes between single jet and dijet. The 
#jet-dependcncc of the G_-modcs for 9j et = {0, 7r/2} is 
displayed in FigJTD] It is seen that the dispersion rela- 



FIG. 9: (Color online) The dispersion relation for the stable 
G--mode for an anisotropy plasma and a jet for 9„ — n/4, 
djet = 0, Vjet = 0.7 and rj — 0.2. The left(right) panel corre- 
sponds to anisotropy parameter £ = 1(10). 
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FIG. 10: (Color online) The dispersion relation for the stable 
G-(G)-mode for the composed of anisotropic plasma with 
single jet and dijet for 8j et = {0, tt/2}, £ = 10, Vj et = 0.7, 
9 n = tt/4 and 77 = 0.2. 



tion for G_-mode is more sensitive to 



'jel 



which is not 



the case for A and G+-modes, The dispersion relation 
for G_-modcs for various 8 n is plotted in FigfTTl We 
find significant dependencies on 9 n . The collective mode 
for G_-mode is also marginally sensitive to the jet veloc- 
ity(see FigfE 



B. Unstable modes 

In the static limit, for small 9 n some of the scales ap- 
pearing in the dispersion relations are negative implying 
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FIG. 11: (Color online) Same as Fig[lO]for £ = 10, « je t = 0.7, 
Ojet = 0, 77 = 0.2 and different 0„ = {tt/4, tt/2}. 
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FIG. 13: (Color online) The growth rate F of the unstable 
mode of single jet and dijet for two different anisotropy pa- 
rameter £ = {1, 10} and 6j et = 0, r] = 0.2 and Vj et = 0.7. 
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FIG. 12: (Color online) Same as Fig[T0]for £ = 10, 6 n = tt/4, 
Ojet = 7r/3, rj = 0.2 and different Vj e t = {0.7, 0.5}. 



that the whole system is unstable with respect to mag- 
netic instability |T7l . [l8j . It has been demonstrated that 
in case of k || n(9 n = 0), the growth rate of the fila- 
mentation instabilities is the largest [T3, El, EH ■ In such 
case 7 and n 2 vanish identically and the dispersion rela- 
tions (henceforth called a and (5 modes) for the unstable 
modes simplify to 

lu 2 -k 2 - a(u) - a'(uj) - a"{uj) = (27) 
w 2 -/3(w)-/3'(w)-/3"(o;) = (28) 



In the numerical solutions of the above equations we have 
checked that unstable mode exits only for a mode both 
for single jet Q and dijet. Let us now concentrate on the 
dispersion relations for the unstable a mode for which we 

shall, first assume that 9j et = 0. It is numerically checked delineated in FigQ3] both for single jet and dijet. For 



FIG. 14: (Color online) Same as FiglTHIwith £ = 10, v jet 
0.7, n = 0.2 and 8 jet = {tt/4, tt/2}. 



that u) is purely imaginary, i.e. lu — iT with T real valued. 
We first show our results for two values of the anisotropy 
parameter £ = {1, 10} and for fixed values of the other 
parameters in the model. This is depicted in Fig|l3] The 
growth of the instability in case of dijet is more compared 
to the single jet case for a fixed value of the anisotropy 
parameter. However, like the case of single jet results Q, 
the instability first increases and then becomes damped. 
It is also seen that the results are quite sensitive to £. 

For 9 j et ^ the roots of the dispersion relations are 
of the form of uj = a + iT, i. e. these are unstable 
propagating modes. This has been checked numerically. 
The results for the unstable modes with non-zero 9j et is 
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0jct = 7r/4 the results are similar to the case of 9j Ct = 0. 
However, for 6 ]ct = 7r/2, the instability always grows for 
dijet case as well as in the case of single jet which has 
already been reported Q. In this case the modes never 
become damped. It is also worthwhile to note that the 
dijet induced instability for 8 n = ir/2 is more compared 
to the single jet case. 

IV. SUMMARY 

The dispersion relations of both stable and unstable 
collective modes in AQGP due to the dijet propagation 
have been obtained and solved numerically. We then 
compare it with that obtained in case of single jet and no 
jet cases. Depending upon the values of different param- 
eters and for the dijet case considered here, the stable 



modes differ significantly from that obtained with single 
jet. We have shown the dependencies of the dispersion 
relations (both for single and dijet) on various parame- 
ters explained in the text. It has been observed that the 
dispersion relation for stable modes are very sensitive to 
the anisotropy parameter, the anisotropy direction and 
the jet direction with respect to the wave vector. For 
the unstable modes it is found that the growth rate with 
dijet is more compared to the single jet case. In partic- 
ular, for 9j et = 7r/2, the growth rate for the dijet case 
is substantially larger and it always grows like the single 
jet case Q. Thus the introduction of dijet might lead to 
faster isotropization as compared to the single jet case. 
For the special case considered here (k||n) no unstable 
modes exits for /3-modc even with the dijet. However, 
more general case might lead to /3-unstablc modes which 
is worth investigating. 



[1] C. Manuel and S. Mrowczynski, Phys. Rev. D 74, 105003 

(2006) . 

[2] M. Mannarelli and C. Manuel, Phys. Rev. D 76, 094007 

(2007) . 

[3] M. Mannarelli and C. Manuel, Phys. Rev. D 77, 054018 

(2008) . 

[4] M. Mandal and P. Roy, arXiv:1204.2976 

[5] M. Mannarelli, C. Manuel, S. Gonzalez-Solis and 

M. Strickland, Phys. Rev. D 81, 074036 (2010). 
[6] S. Mrowczynski, Phys. Rev. C 49, 2191 (1994). 
[7] P. Arnold, J. Lenaghan, G. D. Moore and L. G. Yaffe 

Phys. Rev. Lett. 94, 072302 (2005) 
[8] A. Rebhan, P. Romatschke and M. Strickland Phys. Rev. 

Lett. 94 102303 (2005) 
[9] S. Mrowczynski, Acta Phys.Polon. B 37 427 (2006). 
[10] P. Romatschke and M. Strickland Phys. Rev. D 68, 

036004 (2003). 



[11] P. Romatschke and M. Strickland Phys. Rev. D 70, 

116006 (2004). 
[12] H. T. Elze and U. Heinz Phys. Rep. 183, 81 (1989). 
[13] E. Weibel Phys. Rev. Lett. 2, 83 (1953). 
[14] P. Arnold, D. T. Son and L. G. Yaffe Phys.Rev.D 59 

105020 (1999) 

[15] B. Schenke, M. Strickland, C. Greiner and M. H. Thoma 

Phys. Rev. D 73, 125004 (2006). 
[16] A. Dumitru, Y. Guo, and M. Strickland, Phys. Lett. 

B662, 37 (2008). 
[17] M. Strickland Brazilian Journal of Physics, vol. 37, no. 

2C (2007). 

[18] P. Arnold, J. Lenaghan and G. D. Moor JHEP 08, 002 
(2003); P. Arnold, G. D. Moor and L. G. Yaffe, JHEP 
01, 030 (2003). 



